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Electronic correlations in the ground state of an idealized 
infinite-layer high-T c compound are computed using the ab- 
initio method of local ansatz. Comparisons are made with the 
local-density approximation (LDA) results, and the correla- 
tion functions are analyzed in detail. These correlation func- 
tions are used to determine the effective atomic-interaction 
parameters for model Hamiltonians. On the resulting model, 
doping dependencies of the relevant correlations are investi- 
gated. Aside from the expected strong atomic correlations, 
particular spin correlations arise. The dominating contribu- 
tion is a strong nearest neighbor correlation that is Stoner- 
enhanced due to the closeness of the ground state to the mag- 
netic phase. This feature depends moderately on doping, and 
is absent in a single-band Hubbard model. Our calculated 
spin correlation function is in good qualitative agreement with 
that determined from the neutron scattering experiments for 
a metal. 
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I. INTRODUCTION 

The theoretical understanding of the microscipic elec- 
tronic properties of the high-T c compounds is still in- 
complete. The only ab-initio methods that so far have 
been applied to these compounds are based on the local- 
density approximation (LDA) within the framework of 
the density functional formalism These fail to de- 

scribe some of the basic properties like the magnetic tran- 
sition or the magnetic correlations (for a review see Q). 
Consequently, simplified models have been used that are 
mostly restricted to a single band of strongly correlated 
electrons, and show a Mott-Hubbard localization transi- 
tion at half filling. These seem to explain some of the 
magnetic properties but their microscopic connection to 
the full Hamiltonian has not yet been fully established 
(for a review see 

Here, we present the first application of the local 
ansatz (LA) to these materials. The LA is an ab- 
initio method for the treatment of the correlated elec- 
tronic ground states of solids |]||- It contains no 
homogeneous-electron-gas like approximation wherever, 
and consequently has no problems in overlooking mag- 
netism. The LA yields not only ground-state energies 
or densities but also detailed correlation functions. In 
particular we present the detailed intraplanar correla- 
tion features relevant for all high-T c compounds. Of spe- 
cific interest are the magnetic correlations. We compare 
the frequency integrated momentum dependent inelas- 



tic magnetic neutron scattering intensity measured for 
LaQ.ssSro.isC^Oi [Q to the LA results and connected 
it to a specific correlation. 

The LA is similar to Quantum Chemistry (QC) meth- 
ods, which provide a satisfying description of electrons 
in small molecules. It allows to extend the QC-accuracy 
to solid calculations. Like most of these methods, the 
LA adds correlations as corrections to a single-particle 
self-consistent-field (SCF) ground state obtained from a 
Hartree-Fock (HF) calculation where the electrons are 
described in a restricted single particle basis of Gaussian 
type orbitals (GTO's). The HF-computation for the solid 
is performed by the program Crystal92 ||. At present, 
the LA is the only ab-initio correlation scheme available 
that makes use of this HF program.. Unlike the QC- 
methods, the LA allows a quantitative treatment of elec- 
tronic correlations for solids, independent of their nature, 
i.e. whether they are insulators or metals. The LA can 
do so because it does not attempt to cover the complete 
spaces of one- or two-particle excitation operators as QC- 
methods usually do, but considers the local character of 
the relevant correlations from the very outset. It can 
be seen as an apropriate generalization of the Jastrow- 
ansatz |^] to inhomogeneous systems. 

Every correlation operator in the LA scheme has a very 
specific meaning. It is constructed from pairs of local 
orbitals, each of which is connected with a single atom. 
Consequently, all the incorporated correlation corrections 
are separated into those on single atoms and those inbe- 
tween atom pairs. The full correlation treatment can be 
segmented and partitioned with the individual atoms as 
the smallest available subunit. The correlation opera- 
tors and their treatment are essentially independent of 
the nature of the SCF-ground state, i.e. whether this 
be metallic or otherwise. Such a restrictive choice for 
correlation operators leads to a strong reduction in the 
correlation-operator space and thus substantially facili- 
tates computations. Necessarily, it admits a small loss of 
the correlation energy available in a full treatment within 
a given basis set. From previous calculations, this loss is 
known to be only 2%, independent of the system size 0. 

The LA was used before for extended molecules like 
C*6o [0 , three dimensional semiconductors [ll| and ionic 
insulators fl^l , as well as one-dimensional (polyacety- 
lene |l3|]), two-dimensional (graphite (l4|), and three- 
dimensional metals (Li (l^]). The calculations presented 
in this work concern the first application of the LA to a 
metallic transition-metal compound. 
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From the ab-initio results, also model Hamiltonians 
can be derived that are based on atomic degrees of free- 
dom. The ab-initio correlation functions obtained from 
the LA allow in particular to unequivocally determine the 
model interaction parameters |T^ |. This feature enables 
us to extend correlation calculations to problems that are 
still out of reach for ab-initio LA calculations, and fur- 
ther to make comparisons to models determined by other 
methods. This is particular relevant to the high-T c com- 
pounds, whose properties are often adressed by means of 
model calculations. 

Models based on atomic degrees of freedom consist of 
selected sets of orthogonalizcd atomic orbitals i that are 
related by hopping terms. The interaction part of the 
Hamiltonian, H™° d , is usually restricted to local interac- 
tions 
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for two electrons in the same orbital i. Each interaction- 
energy parameter Ui dominantly influences a particular 
correlation function of the correlated ground state ^ cor r 
for such a model, namely C orr\n-i-\riii\^ corr ) , or equiv- 
alently the change of this correlation function due to cor- 
relations, 



Ai(corr) = (* corr KT n i|.l*c rr) - (^SCfKt^iI* 



SCFI 

(2) 

Here, ^scf represents the SCF-ground state of the 
model. As for the ab-initio treatment, the model ground 
state is also computed by means of the LA. Ai(corr) 
equals zero for Ui = and rises continuously with Ui. 
When multiplied with Ui, it is a measure of the interac- 
tion energy. This is the relevant relative quantity, when 
symmetry is broken due to the atomic interactions, and 
a magnetic or structural phase transition occurs. In such 
a case, states are compared that differ in those atomic 
fluctuations. These energy costs due to local charge fluc- 
tuations are also relevant in the context of Compton scat- 
tering, secondary or shake up peaks in photo emission or 
core spectroscopy. For the transition metals, it turns out 
that the same model interactions are needed for the de- 
scription of all these properties p^-[20|] . 

With the unequivocally defined orthogonalized atomic 
orbitals available in the ab-initio calculation (see section 
H), this same correlation function can be determined 
from the ab-initio calculation. The model interaction Ui 
can thus be fixed by demanding that the correspond- 
ing model correlation corrections Ai(corr) agree to the 
same ab-initio quantities. This connection can also be 
used to analyse screening details entering such a model 
interaction. By adding stepwise particular screening con- 
tributions in the ab-initio calculations, we will deter- 
mine how the model interaction arises, starting from the 
bare Coulomb interaction. Such an investigation was in 
the past performed for 7r-electron interactions in organic 
compounds Jig]. 



For the high-Tc materials, ab-initio correlation func- 
tions have not yet been available. Model interactions, 
however have been computed from the LDA calculations 
in a different approximation, by freezing specific charges 
on individual atoms and relaxing the environment. For 
the most extended model that has been used in this con- 
text, a three band model, the two interaction param- 
eters, namely an effective local interaction Ud between 
the Cu3<i-electrons, and an interaction U p between the 
02p-electrons have been obtained in this way [pl|-p3|. 
These interactions have then been used to extend the 
LDA calculations to magnetic properties (by means of 
the so-called LDA+U method or by the self interaction 
correction (SIC) calculations) (^4). Model calculations 
based on these interactions yield photoemission results 
that are in very good agreement with experiment p5[ . 
This agreement is in contrast to the case of in partic- 
ular the middle of the transition metal series where ef- 
fective interactions obtained by the LDA do not match 
the interactions needed for the description of the above- 
mentioned properties of these systems (for a comparison, 
see 16 ). For the high-T c compound that we shall deal 
with, we will make a comparison between the effective 
interactions obtained by the LA and by the LDA. It is 
the first of its kind because for the systems treated so far 
by the LA no LDA interactions are available. The com- 
parison will shed additional light on the transition-metal 
case. 

For simpler single-band models, interactions have so 
far been mostly guessed, or obtained by means of fits to 
specific experiments. Often, they are deduced from the 
LDA three-band models. However, usually the strong in- 
teraction found by the LDA is the only transferred quan- 
tity. The Cu3<i a .2_j ( 2-occupation of 1.5 of the LDA |24] 
was usually replaced by occupations smaller than 1.2 W. 
Only under such conditions, a Mott- Hubbard scenario 
applies iffij . 

This paper is organized as follows. In section 0, a short 
description of the LA and of its possible shortcomings in 
connection with its application to high-T c compounds is 
given. Section III contains a detailed discussion of the 



ab-initio computation and an analysis of the correlation 
functions. These calculations are performed for an ide- 
alized high-T c compound, the so called infinite-layer sys- 
tem. Whereever possible, comparison to the LDA results 
is made. Sections |y| and |y| deal with the computation 
of the model interactions. The former section contains 
in particular a detailed account of the screening mecha- 
nisms for the 3d transition-metal interaction. In section 
VI, the dependencies of the most important correlation 
features on band filling are investigated. This analysis is 
made on the model level. The ab-initio program Crys- 
tal92 that is used for the SCF calculation is restricted 
to integer electron occupations per unit cell, and a small 
change of band filling on the ab-initio level would de- 
mand large unit cells. Finally in section VH, a compari- 



son to neutron-scattering results for a metallic compound 
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is made, and the experimentally found inelastic scatter- 
ing is connected to a particular correlation obtained by 
the LA calculation. A su mmar y of our work is presented 
in the concluding section 
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II. COMPUTATIONAL DETAILS 

The aim of the paper is a quantitative understand- 
ing of the electronic correlations relevant to all high-T c 
compounds, namely those in the individual CuO-plane. 
To simplify the computations, the idealized compound 
SrCu02 is treated. It is a so called infinite layer com- 
pound with CuO-plancs separated by layers of Sr ions. 
In the planes, the Cw-atoms form a quadratic net with 
a separation of 3.925 A, with the O-atoms at equal dis- 
tances inbetween neighbor Cu-atoms. The stacking dis- 
tance, and the perpendicular Cu-separation amounts to 
3.43 A. The Sr -ions have equal distance to four Cu- 
atoms in the two neighbor planes each. This compound 
has the smallest possible unit cell containing four atoms. 
The uppermost valence band is half filled. Therefore, 
the compound is expected to order antiferromagneti- 
cally. Our interest , however, is in the correlations in the 
metallic state. To represent a doped metallic compound 
would require a very much larger unit cell. Instead, the 
metastable non magnetic state for the small unit cell is 
used as a starting point. This approximation is justi- 
fied because it is known from LDA calculations that the 
energy bands and Fermi surface of the metastable half- 
filled metallic state are very close to the ones for the true 
metallic compounds |p|,p7|. 

In a first approximation, this metastable metallic state 
is obtained from a restricted, non symmetry broken 
Hartree-Fock calculation for this compound, performed 
with the ab-initio program Crystal92 ||. For the Cu 
and the O-atoms, good all-electron-GTO-basis sets are 
used. For Cu, this is a modified (14,11,6) Ahlrichs ba- 
sis |f28| , contracted to (6,4,2) orbitals. From the original 
basis, the outermost diffuse functions were removed, and 
the next exponents adjusted and reoptimized. For O. 
this is a (11,7) Huzinaga basis p9| , whose outermost ex- 
ponents were contracted as was done before jn],[l2| > phis 
a set of d-orbitals. While the basis sets for the atoms 
in the planes are of good quality and promise results for 
the valence electrons close to the Hartree Fock limit, this 
is not the case for the Sr atoms. The latter are rep- 
resented by a large core pseudopotential and a single 5s 
orbital p3]. Here, also the outermost diffuse basis orbital 
was removed. There is no need for such a treatment, but 
due to this choice, the charge distribution and correla- 
tion analysis can be definitely restricted to the degrees of 
freedom within the plane. 

From the SCF calculation, the metallic single-particle 
ground state |\&scf) is obtained. Its Fermi surface is 
identical to the one of a LDA-calculation for the same 
compound p7[. Also the uppermost energy band is sim- 



ilar to the equilavent LDA band, except of an additional 
homogeneous spreading by almost a factor of two due to 
the non-local and non-screened exchange. The non-local 
exchange also causes the lower lying bands to be more 
separated from the uppermost half-filled band than in 
the LDA calculation. A presentation of such HF bands, 
obtained in a somewhat different basis, is found in Ref. 
0- 

In a next step, correlations are added by the LA. Here, 
the following variational ansatz is made for the correlated 
ground state: 



s = 
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I^scf) 
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The ry's serve as variational parameters. The rii a and Sj 
are density and spin operators for an electron in the local 
state aL , represented by the orbital 
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where the fi(r) are the (GTO like) basis orbitals. The 
operators have an obvious meaning. The first operator 
riiiriii, for example, when applied to I^scf), projects out 
all configurations with two electrons in orbital gi(r). In 
connection with the variational parameter r\ v , as in eq. |J, 
it partially suppresses those configurations. Similarly, 
the operators mrij describe density correlations between 
electrons in local orbitals gi(r) and gj(f). For the homo- 
geneous electron gas, an ansatz with these two kinds of 
operators leads to the Jastrow function 0. The opera- 
tors s*j • Sj generate spin correlations. The fourth kind 
of operators is of the form of [O v , Hq] _ , where Ho repre- 
sents the single-particle Hamiltonian. In comparison to 
the first three kinds of operators which look like partic- 
ular interaction contributions, these operators refine the 
ansatz with respect to the band energy of the electrons 
involved Within the computation, the original oper- 
ators of eq. P are modified by subtracting the contracted 
contributions in each of them. The corrected operators 
when applied to |v&scf) contain only two-particle excita- 
tions, and the corrected last kind of operators in eq. |5| 
covers local single particle excitations, i.e. it allows for 
changes in occupations. 

The variational parameters r/^ are chosen to optimize 
the energy 
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In the last equation, the subscript c indicates that only 
connected diagram contributions are summed up. This 
expression cannot be evaluated exactly. The standard 
approximation is an expansion in powers of 77, up to sec- 
ond order, 



Eg — Escf + E COII 



(9) 
(10) 

(11) 



Here, (A) means the expectation value of the operator 
A within |\&scf)- This approximation works only if the 
correlations are sufficiently weak. Disregarding the re- 
duced subspace of correlation operators, the approxima- 
tion used so far corresponds to a Linearized Coupled 
Cluster Singlet and Doublet (LCCSD) treatment @. It 
can be extended to a CCSD treatment 1 33 1 . 

The local orbitals in eq. ^ are connected to a single 
atom only and are built from its basis orbitals. This is 
the essential approximation of the LA. 

In the present application, only atomic orbitals are 
constructed. These are uniquely determined from the 
SCF-ground state by the condition that they are built 
from basis orbitals on the respective atoms only and that 
they contain a maximal part of the occupied space. The 
resulting orbitals are next Lowdin-orthogonalized to each 
other. More localized subatomic orbitals that are usu- 
ally included in applications of the LA are not used in 
this first application for a metallic high- Tc compound. 
Therefore, only interatomic correlations are treated, i.e. 
correlations that are described by operators built from 
the atomic orbitals. Shorter range or intra-atomic corre- 
lations, as well as particular polarization correlations are 
not covered. From previous experience, it is known that 
such contributions are not very relevant for the topics of 
interest here. Estimated corrections due to the omitted 
correlations will be given were they are nonnegligible. 

In the LA, the correlations are taken into account in- 
crementally. The correlation energy is exactly expressed 
as an incremental sum over contributions from different 
sets of atom clusters 



E r , 
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m=l [ ji j 2 jm ) 

with jx ^ j 2 ^ . . . ^ jm (12; 

where the Aj n denote atoms, on and between which cor- 
relation operators are formed and it holds , for example, 

(E COII (A 1 A 2 )) i = E COII (A 1 A 2 ) - S corr (Ai) - E COII (A 2 ) 

(13) 

i.e. the increments ()j include only the changes of the 
correlation energy due to the extended set represented. 



Translation invariance and the particular local symmetry 
are easily included by performing the above summation 
only over the subset of symmetry inequivalent clusters. 
For every computation, the exact solid single-particle ex- 
pectation values are taken. 

For the coverage of the interaction part, the local na- 
ture of the correlation operators allows a drastically sim- 
plifying reduction. In a finite basis set per atom repre- 
sentation, the interaction is represented by a fourth order 
tensor of basis interaction matrix elements whose indices 
extend over the involved basis orbitals. The generation 
of this tensor and its handling are the limiting steps in 
a correlation calculation. For the particular correlations 
on the set of atoms treated, the required tensor can, to 
a very good approximation, be restricted to these atoms 
plus all their nearest neighbours. This restriction makes 
all required computations easily feasible. Possible correc- 
tions due to lacking matrix elements are included in com- 
putations extending over larger clusters. The two largest 
clusters for which explicit ab-initio correlation calcula- 
tions were performed are depicted in Fig. |[ The larger 
one consists of five active Cu atoms and four active O 
atoms. The basis interaction matrix elements are com- 
puted for the whole Cu^Oiq cluster. As will be demon- 
strated, the most relevant informations can be satisfacto- 
rily obtained from correlation calculations extending up 
to this size of clusters. 

Two approximations made for the handling of the LA 
need further discussion in connection with the applica- 
tion to a high-Tc compound. The one is the restric- 
tion to weak correlations while the high-Tc compounds 
are usually connected with a Mott-Hubbard transition. 
The LCCSD approximation, in which the LA is com- 
puted, fails for the strongly correlated half-filled band 
case. However it does so in a controllable way. For too 
strong correlations, the correlation corrections turn too 
large and lead to negative density correlations functions. 
The criterion of positive density correlations can be taken 
as an indication whether the LA results are still mean- 
ingful. Away from half-filling, the LA behaves better, 
and for an almost empty band, it applies even for diverg- 
ing interaction. Also, a more quantitative test can be 
made. Individual correlation corrections like the one due 
to a single correlation operator n^riii can be computed 
variationally, i.e. exactly. Here, a comparison with the 
result of the eqs. |ll| in the same operator subspace can 
be made, and the overestimation of the correlation ex- 
pansion can be quantified. Such a variational calculation 
restricting to two-particle excitations has no meaning for 
the treatment of the full extended system (N) due to 
lack of size consistency. The resulting correlation energy 
would scale like \/N. When such a variational compu- 
tation is extended to more then a single operator, it can 
only give a lower limit to the correlation results. Ear- 
lier ab-initio calculations with the LA for finite Cu — O 
clusters |34|] and experience with succesful LA mod el cal- 
culations for the transition metals themselves ItS 
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have already indicated that the range of applicability of 
the LA might well extend to the high-T c materials. 

The second approximation made is the restricted SCF- 
ground state that is used as a starting point, although the 
particular system chosen is known to be antiferromagnet- 
ically ordered. The LCCSD approximation is sensitive to 
antifcrromagnetic order. This is in contrast to standard 
perturbation expansions, and can be seen by resolving 
the LCCSD equations (pi]), leading to 

E COII = -J2{OtH)({O^HO) c )-^{OlH) . (14) 

The denominator contains the exact two-particle excita- 
tion energies. If the restricted SCF-ground state turns 
unstable, its susceptibility diverges for a particular wave 
vector, and consequently two-particle excited states must 
exist with energies degenerate to or even lower than the 
SCF ground state energy. With sufficiently extended sets 
of correlation operators, the matrix {OHO) is no more 
positive definite, and the scheme turns instable. From 
such a calculation, also the smallest set of correlation op- 
erators may be determined that leads to instability. In 
particular, a lower limit for the size of stable magnetically 
ordered domains can be obtained. This holds true as long 
as the phase transition from the metal to the insulator 
as a function of occupation is second order. If the transi- 
tion is first order, then the computation in the metastable 
metallic state is still relevant for the doped metallic state 
but might lack informations about the magnetic order. 



There are two different kind of relevant longer range 
contributions. One arises from correlation operators be- 
tween neighbor Cu and O atoms. Here, no specific contri- 
bution is dominating. Rather the very local atomic corre- 
lation hole generated by the atomic operators is smoothly 
extended, adding ten percent of the on-site correlation 
energies. The second kind is connected with spin cor- 
relations between different Cu-atoms, and in particular 
with those inbetween electrons in the C , M3d :E 2_j ( 2-orbitals. 
These dominate the nearest neighbor Cu — Cu correla- 
tions and are exclusively responsible for the longer range 
terms. The neighbor Cu — Cu correlations will be later 
elaborated in more detail. The longer range correlations 
are connected with the eventual formation of long range 
magnetic order. Next nearest Cu — Cu contributions to 
the energy are as large as nearest neighbor contributions, 
indicating that here problems with the metastability of 
the nonmagnetic SCF-ground state begin to show up. 

While all shorter range correlations were fully or al- 
most completely converged with respect to the series of 
clusters selected, this does not hold true anymore for the 
third nearest neighbor Cu — Cu-contributions. However, 
divergence doesn't yet appear. This indicates that sta- 
ble antiferromagnetic correlations at half-filling need to 
coherently extend over domains larger than the clusters 
selected for the correlation computation. 

When added, all these correlations represent an energy 
gain of roughly heV per unit cell. All these correlations 
are due to binding, and the resulting correlation energy 
amounts to a large fraction of the total binding energy. 



III. RESULTS OF THE CORRELATION 
CALCULATIONS 

We will next present the ab-initio results of the LA 
in separate chapters for correlation energies, charge dis- 
tributions and particular correlation functions. In all 
cases, partial correlations are consecutively added, start- 
ing from atomic terms, and extending up to the longest- 
range correlations included, namely the ones between 
third-nearest-neighbor Cu-atoms. 

A. Detailed correlation energies 

In Table |, the contributions of the different classes of 
operators to the correlation energy are displayed. The 
sets of operators are grouped into those on individual 
atoms and those inbetween different atoms. For the lat- 
ter cases, also the fraction of spin operator contributions 
is given. As expected, the largest overall energy gain 
is due to the on site or atomic correlations. Here, the 
largest part is from the Cu — 3d x 2_ y 2 operators. How- 
ever, almost twenty percent of the on site contributions 
is connected with a charge transfer that will be discussed 
in more detail later. 



B. Partial charge distributions 

The 

partial charge distributions n^) = (*| J2 a riia]^) are 
presented in table || for different states ^. The first 
row contains the values for W = 't'scF- When added, 
the partial occupations reach the number of valence elec- 
trons up to 0.02. This indicates the good quality of the 
computed orthogonal atomic orbitals. The occupation of 
the Cit3d-orbitals is very close to the estimate for the 
solid, obtained from earlier finite cluster HF-calculations 
, except for the Cu4s, 4p-occupations. In the earlier 
calculation, these came out smaller for two reasons. The 
one is that a basis set was used that lacked the most ex- 
tended exponents used here for the 4s,p-orbitals. This 
restriction was made to avoid artifacts, resulting from the 
large negative charging of the small clusters treated. The 
second reason is that in the earlier calculation, the 4s and 
4p-orbitals were Schmidt orthogonalized to the 02s, 2p- 
orbitals, while here all orbitals are equally treated by a 
mutual Lowdin orthogonalization. 

With the addition of correlations, a relatively large 
charge transfer occurs. Ultimately, it is a charge trans- 
fer mostly from the Cu3d x 2_ y 2 -orbitals into the 02p- 
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orbitals. However, for its understanding it is neces- 
sary to progress stepwise. A first step is the addition 
of atomic correlations which lead to a large correlation 
energy gain. The dominant charge transfer due to the 
atomic correlations is from the Cu3c? 2 ,2_ i) 2-orbitals to the 
CuAs, 4p-orbitals, followed by a secondary redistribution 
from the Cu4s, 4p-orbitals to the 02p-orbitals. Over all, 
0.18 electrons are removed from the Cu3dj.2_ a 2-orbitals, 
and put into the Ctt4s,p-shell (0.13) and the 02s,p-shells 
(2x0.03). This charge transfer was not detected in the 
earlier cluster calculation for the poor Cit4s, p-basis pi| . 
More than 80 percent of this charge transfer arise from 
the inclusion of the operators n^riii for the Cu3d x 2_ y 2- 
orbitals, the remaining part stems from the same opera- 
tors for the 4s,p-orbitals. 

This charge transfer due to atomic correlations is 
closely related to the negative magneto-volume effect 
known from transition metals as will be shown next. 
The correlation induced charge transfer detected here 
partially corrects an inverse exchange induced charge 
transfer. The dominant exchange contribution of rele- 
vance in this context is from atomic interactions U a t(i) 
and is written as E at (exch) = -Y,i U at{i)n 2 ia {^ hf) 
where the sum runs over the i atomic orbitals in the 
unit cell, and n i(J (^HF) indicates the HF-occupation of 
orbital i per spin. To simplify the discussion, we re- 
strict to two different atomic orbitals, and assume that 
they have the same interaction U but an occupation n\ a 
very much larger than an occupation ri2 a . A charge 
transfer z from ni into n\ leads to the exchange en- 
ergy gain AE(exch) = —2zU(ni a — nia)- Therefore, 
the HF-exchange enhances differences in charge distri- 
butions. For the case treated here, the HF-exchange is 
responsible for a charge transfer from the little filled 4s, p- 
orbitals into the strongly filled 3<i-orbitals. Also a trend 
for a similar charge transfer from the 4s,p-orbitals into 
the strongly filled 02s,p-orbitals is expected, but should 
be restricted by Hartree (or essentially Madelung) con- 
tributions. The latter terms don't influence intra-atomic 
charge transfer. No charge transfer between the Cu3d 
and the 02p-orbitals is expected from atomic exchange 
as long as these orbitals have a similar occupation and 
not very different atomic interactions. When on-site or 
atomic correlations are included, a sizable fraction of the 
exchange induced charge transfer is undone. This is the 
origin of the charge redistribution found in the present 
computation. Another way to undo the exchange induced 
charge transfer is to turn the system magnetic. This cor- 
responds to a maximal atomic correlation. Consequently, 
a charge transfer must come into play when magnetism 
in systems with very differently filled subshells is con- 
cerned, as for the case of itinerant ferromagnetism of the 
3d-metals. For Ni, for example one would expect a some- 
what larger filling of the 3<i-orbitals for the non magnetic 
state than for the ferromagnetic state at the cost of the 
4s, j>occupation, because contrary to the fully magnetic 
case, the electrons are not completely correlated in the 
non magnetic state, as is well known pq. 



Such a charge transfer becomes relevant for the 
magneto- volume effect. A theoretical description of mag- 
netism that is restricted to a particular shell (like the 
set of 3d orbitals) implies that the antibonding orbitals 
are more populated in the magnetic case no matter how 
strong correlations are. In such a case, the volume al- 
ways increases with magnetic order. When more than 
one shell is involved, the just mentioned charge transfer 
comes into play. For Ni, this is from the bonding 4s, p- 
orbitals into the (few and antibonding) empty orbitals 
of the 3d-shell. It implies an negative magneto volume 
contribution. Ni actually displays a total negative mag- 
netovolume effect |3(| . The latter can not be understood 
within a description restricted to the 3<i-orbitals Jl8| , but 
can only be explained by a not fully screened exchange 
induced charge transfer. On the HF level, this exchange 
induced 3d — 4s-charge transfer was first proposed many 
years ago as the origin of the negative magneto-volume 
effect fl37| . 

It would be of interest to find out whether also in the 
case of the high-T c compounds a negative magnetovolume 
effect exists. With the new version of the Crystal pro- 
gram p8fl , which allows for unrestricted HF-calculations, 
such an investigation will become feasible from the the- 
oretical side. 

When neighbor correlations are included, then an ad- 
ditional charge transfer of the same magnitude as the 
one due to on-site correlations occurs. It is dominantly 
from the Cu3d,,2_ a 2-orbitals to the 02pb-orbitals, and 
is due to a particular spin correlation between neighbor 
Cu-sites that will be discussed later. The longer range 
contributions that were covered by the present compu- 
tations lead to a further but small transfer of the same 
kind. A similar but somewhat smaller charge transfer 
had been found in the earlier finite cluster calculations 
. It should be noted that this secondary charge trans- 
fer is connected with a small correlation energy gain. 

The exact occupation of the Cu3d x 2_ y 2 orbitals is 
of very much interest. It is directly connected to the 
measured moment of the magnetic ground state and 
plays a crucial role for models used for the high-T c com- 
pounds. These usually assume that the occupation of 
the Cu3d x 2_ y 2 -orbital is not very different from 1.0. The 
present calculations indicate that in fact different corre- 
lation mechanisms bring the occupation into this range. 
To validify these findings, a short discussion on the pos- 
sible deficiencies of the presented computation shall be 
given next. 

While the result of the HF-calculation can be assumed 
to be close to the HF-limit, this does not hold true for the 
correlation treatment. A first possible error is connected 
with the weak correlation approximation. From the ex- 
act treatment of a particular on-site correlation plus the 
connected charge transfer, it can be estimated that the 
on-site correlation correction of the Cu3d x 2_ y 2- occupa- 
tion is overestimated by 10-15 percent. For the charge 
transfer arising from longer range correlations, no error 
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estimate can be made. Two further corrections are ex- 
pected. The one is the influence of shorter than atomic 
range correlations, the other is the screening particularly 
of the 3c£-orbitals due to the Cu core - and here especially 
due to the 3s and 3p orbitals. Since the calculation pre- 
sented here is the first ab-initio correlation calculation 
ever for a metallic transition-metal compound, no refer- 
ence results exist. Even comparable detailed calculations 
for small clusters are lacking. From atomic calculations, 
some estimate for corrections can be gained. Such cor- 
rections are expressed in terms of energy differences. In 
HF-approximation, the excited 3dg4s2 state of the Cu 
atom is 0.4 eV higher than the ground state, while ex- 
perimentally the difference amounts to 1.5 eV A 
similar correction is found when the ground and excited 
states of the TVi-atom are concerned 1 40 . This correc- 



tion is exclusively due to short range correlations, core 
polarization effects, and relativistic corrections, none of 
which were included in the presented calculation. The 
influence of these corrections to the presented results can 
be estimated by lowering the diagonal atomic 3d-energy 
level by 1 eV. For the model discussed below, this leads to 
a charge transfer into the Cu3(i ; j,2_j ( 2-orbitals of roughly 
0.05. The lacking short range correlations alone would 
have a somewhat larger influence than this total shift 
but are counterbalanced by relativistic corrections that 
favor the 4s-electrons Ej]]. Adding these corrections, the 
final estimate for the Cw3(i :I .2_ J( 2-occupation amounts to 
1.22 ± 0.07. For reasons discussed above, it is smaller 
than the value deduced from the earlier cluster calcula- 
tions |34). Future applications will hopefully reduce the 
uncertainty in the present LA results. 

It is of interest to compare this LA result to an LDA- 
charge analysis. Table [0] also contains LDA results that 
were obtained for Y BC12CU3OQ.5 |42|. The latter com- 
pound is not a so called half-filled system, and the re- 
spective charge analysis represents a lower limit to the 
half-filled case. In the referred publication, only inte- 
grated occupations for complete shells were given. Also, 
the underlying charge analysis was performed differently. 
This might lead to sizable deviations when the more 
delocalized orbitals are concerned but is hoped to lead 
to comparable results for the very localized 3<i-orbitals. 
From the global 3d occupation, not much can be con- 
cluded about the Cu&d x i _ y i -occupation. However, an 
LDA-calculation performed for the system treated here 
at half-filling leads to a Cu3<i I 2_ J( 2-occupation of 1.55 
p3| , and other high-T c compounds at half-filling are usu- 
ally mapped by Cu3d ;I ;2_ y 2-occupations of 1.5 |25| . These 
LDA values are very close to the HF-results but differ 
from the correlation result and from the final estimate. 

This deviation of the LDA result from the LA occu- 
pation is expected to result from deficiencies of the used 
homogeneous electron gas approximation. It is plausi- 
ble to conclude that the specific neighbor Cu spin cor- 
relations leading to a charge transfer of 0.17 are not at 
all covered by the LDA. Such a simple connection can 



not be made with respect to the charge transfer caused 
by on-site exchange terms and the correlation compensa- 
tions. It is known that the LDA is not able to describe 
anisotropic exchange contributions. In particular, the 
LDA is not able to produce the negative magneto- volume 
effect for Ni . This indicates that it lacks exchange in- 
duced transfer and partial correlation compensation, but 
no large overall error is expected on the atomic scale. A 
very rough error estimate of the charge distribution due 
to LDA deficiencies on the atomic scale can be made us- 
ing an analysis of LDA results for two-atomic clusters. 
These indicate that the 3d-orbitals are too attractive in 
comparison to the 4s-orbitals. Expressed in diagonal en- 
ergies, a correcting shift of 1 eV was computed j|5). Such 
an atomic LDA correction is similar in size to the joint 
correlation/relativistic correction of the HF-energy dif- 
ferences for the atoms but has a different prefactor. It 
leads to a charge transfer of 0.05 out of the Cu3d x 2_ y 2- 
orbital. 

Over all, the LDA seems to overestimate the occupa- 
tion of the Cu3d x 2_ y 2-orhitals by roughly 0.3, a large 
fraction of which is explained. It should be kept in mind, 
however, that the atomic orbitals are differently defined 
in both methods. This causes some uncertainty. In the 
future, it will be possible, to perform LDA calculations 
with a new version of the Crystal program [ j38f , and to 
analyze the results by the LA routines, so that at least 
this last uncertainty can be removed. 

The mutual influence of correlations and charge redis- 
tributions is of relevance to ab-initio methods that try 
to adress correlations with Monte Carlo schemes. Varia- 



tional Monte Carlo calculations |46| as well as diffusion 
Monte Carlo calculations (for a review, see jl7|) rely on 
a good trial state. For the first method, the charge dis- 
tribution of the trial state is usually frozen to avoid the 
optimization of very costly external variational parame- 
ters, while the second method is restricted by the frozen 
nodes of the ground state wave function. The findings 
of the LA calculation indicate that a better trial state 
than the so far always selected LDA ground state wave 
function might be needed. 



C. Atomic correlations 

Next, we discuss the individual atomic correlations and 
their strength. The average occupation of orbital i in the 
correlated ground state is defined as rii — J2 a n io-(^corr)- 
The charge fluctuations within orbital i, An?, are given 
as 



= n i (l-(^-))+2A i (corr) 
= 2(A l (HF) + A t (corr)) 



(15) 
(16) 
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They are separated into the HF-charge fluctuations, 
Ai(HF), and the correlation corrections Ai(corr) (see 
Eq. ||). The former are defined as the charge fluctu- 
ations in a fictitious single-particle state that has the 
charge distribution of the correlated ground state. These 
charge fluctuation represent the electronic mobility. For 
a single-particle state, it holds that 



(17) 



where represents the density matrix elements per spin 
to all other orthogonal atomic orbitals. The kinetic en- 
ergy gain due to derealization of the electrons in this 
state is proportional to ^2rj\Pij, with the summation 
(j) restricted to nearest neighbors of i. Consequently, 
the reduction of Ai(HF) due to correlations, A, (corr), 
gives also a rough measure of how much band energy is 
lost by the correlations. 

Charge fluctuatons can only be completely frozen out 
for half-filling, i.e. for ni — 1.0. In all other other cases, 
there is a maximal correlation reduction, Ai(corr,max), 
which amounts to 



Aiicorr, max) — 



(f) 2 
(1- 



for Hi 
for 7i, 



< 1 
> 1 



(18) 



The relative correlation strength Q is defined as £ 

Ai (corr) 



Ai(corr,rnax) ' 

All these quantities are given in table [II for the final 
result of the correlated ground state. As can be seen, cor- 
relations are strongest for the Cu3c? a; 2_j ( 2-orbital. Never- 
theless, even there, half of the original fluctuations sur- 
vive, indicating that the electrons are still very delocal- 
ized, and that a renormalization of the effective mass due 
to atomic correlations of no more than 30 percent is to be 
expected. Nevertheless, 70 percent of the possible reduc- 
tions are realized. Next in strength are the correlations 
on the 02pb-orbitals. Also here, the correlation strength 
is 0.7 although the reductions amount to only 20 per- 
cent of all fluctuation in this orbital. The correlations in 
all other orbitals are weak. This even holds true for the 
02 S -orbital. 

The correlation strength strongly depends on the in- 
cluded correlations. When restricting to on-site correla- 
tions, the Cu3(i ;1 ,2_j ( 2-occupation is 1.33. Then, it holds 
that Ai(corr) — —0.096, which represents a 85 percent 
reduction. Freezing the Cu3d x 2_ y 2-ch&rge at the HF, 
value, i.e. close to the value of the 02pf,-charge leads to 
a correlation strength of more than 0.90. This will be ex- 
plained later when analyzing these correlation functions 
in the context of model interactions. 

A set of trial variational calculations restricted to in- 
dividual correlation operators was also performed. This 
was done to control the validity of the variational expan- 
sion. When comparing these variational results to the 
variational expansion results, it was found that the corre- 
lations obtained by the expansion calculation, were over- 
estimated by 10-15 percent for the Cit3d.j;2_y2-orbitals, 



but less than 5 percent for all other on-site correlations. 
This small corrections indicate that the expansion is fully 
able to cover these correlations. Consequently, there is 
no evidence that correlations on the atomic scale are too 
strong for a weak correlation expansion treatment. 

There are additional corrections expected from the 
omitted correlations. From earlier calculations for other 
systems, it was found that the longer range correlations 
that were neglected here have no influence on the atomic 
correlation functions (for a detailed explanation, see eg 
p6|). The short range correlations omitted here, how- 
ever, led to a reduction of the atomic correlation correc- 
tions for the 2s,p-orbitals by 10 ± 5 percent Jig ]. Conse- 
quently, they are expected to lead to a somewhat larger 
reduction for the 3d-orbitals which are characterized by 
a somewhat higher average density. When added, a re- 
duced correlation strength of 0.53 ± 0.04 instead of 0.7 
is expected for the C'u3d x 2 _ y 2-orbitals, when an occupa- 
tion of 1.17 is assumed. Such a correlation correction is 
not terribly large but is rather similar to the correlation 
strength obtained for the transition metals from model 
calculations p7[ . 

This correlation strength obtained for the metastable 
haf-filled state should not strongly change when elec- 
trons are removed from the planes. Actually, for such 
a case, the Cu3d x 2_ y 2 -occupation comes closer to 1.0, 
and Ai(corr) is expected to increase while Q decreases. 
Therefore, the results obtained for the half-filled case are 
expected to be representative for relevant dopings. 

Complete magnetic order implies a correlation strength 
of 1.0. Consequently, we would expect a certain addi- 
tional charge transfer 3c?— 4s, p with magnetization. This 
might even lead to a negative magnetovolume effect. 



D. Spin correlations 

Of particular interest are the spin correlation functions 
S(i,j) between the different atomic orbitals , de- 

fined as 



S(i,j) 



(19) 



For the on site terms, it holds in general that S(i, i) = 
l(A l (HF)-A i (corr)). For the SCF- ground state, these 
expectation values don't vanish, and represent the auto- 
correlations of the electrons. They are small except for 
the on site terms and the neighbor Cu — O contribu- 
tions. For the correlated ground state, it turned out that 
only spin correlations inbetween the Cv3>d x i_ y i -orbitals 
on the different atoms are relevant. This is in contrast to 
earlier calculations for small clusters. There, longer range 
spin density wave like correlations occured at which the 
02ph-orbitals participated [pi| . Apparently, these small 
clusters rather represented a 1-dimensional chain than 
the intended plane. 

Table [V displays the correlations S(i,j), where the 
indices («7j) are restricted to Cu3d x 2 _ y 2-orbitSL\s. These 
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correlation functions are computed incrementally. In a 
first step, only on-site correlations are included, and only 
the on-site function is given, then neighbor interactions 
plus the neighbor functions are successively included. 
The enhancement of the on-site terms with inclusion of 
longer range correlations is due to the enlarged charge 
transfer. 

The antiferromagnetic Cw-neighbor spin correlations 
are large and are connected with a sizable charge trans- 
fer from Cu to O that was discussed above. The cou- 
pling of spin correlations to a charge transfer arises be- 
cause spin correlations are maximal for a Cu3d x 2_ y 2 oc- 
cupation of 1. In more detail, the quantity most rele- 
vant for the spin correlation strength is the expecation 
value (siSjH), where i,j represent neighbor Cuid x 2_ y 2 
orbitals. The dominant contributions from the Hamilto- 
nian are the on-site interactions on Cu, U^d, and on O, 
U2 P - When restricting to these interaction, and to the 
largest density matrix elements Pjj , the above expecta- 
tion value is given by 



(sls-H) 



r(l 



Ma)Pi~U 2p Pt) ■ (20) 



Here, rii a represents the SCF-occupation of the Cu- 
orbital per spin, and P\ the neighbor Cu — O density 
matrix element and Pi the neighbor Cu — Cu density 
matrix element. It holds that ni„ > \Pi\ > I-P2I for the 
relevant range of occupations. For the half-filled SCF- 
ground state, the first term of Eq. ^0| is one order of 
magnitude larger than the second. The correlation in- 
duced charge transfer results probably because the ma- 
trix element (O u H) is enhanced by a charge transfer from 
the Cit-orbital to the O-orbital. Due to such a transfer, 
the first part increases, while the second part strongly de- 
creases. The charge transfer stops when the Cuid x 2_ y 2- 
occupation reaches 1. 

The ab-initio results show that the charge transfer due 
to the secondary spin correlation is large. This suggests 
that the ground state must also be relatively instable 
with respect to any other external disturbance that prof- 
its from such a charge transfer. May be this instability 
also contributes to the lattice instability, and in particu- 
lar to the large buckling found for this compound 0] . In 
view of this sensitivity, it is even more astonishing that 
the very much larger on-site correlations did not lead to 
a sizable Cu — O-charge transfer. 

These neighbor correlations are quite strong, stronger 
than needed to counterbalance the change in the wave 
function due to on site correlations. For a singlet state 
f it always holds that X^-W^il*) = °- For the SCF " 
ground state, the on-site terms are counterbalanced in 
part by antiferromagnetic neighbor Cu — O spin correla- 
tions, while the remaining correction is rather long range. 
When correlations up to neighbor Cu atoms are included, 
then the neighbor Cu — Cu correlations alone more than 
counterbalance the correlation enhanced on-site terms. 
This can be viewed as a quite sizable attraction of elec- 
trons of different spin on neighbor sites even in the ab- 



sence of longer range antiferromagnetic order. In fact, it 
will be later demonstrated that these correlations depend 
little on doping. 

While the nearest neighbor correlations represented in 
table IV are converged with respect to the treated cluster 
sizes (as long as no longer range correlation operators 
are added), this does not quite hold true for the second 
nearest neighbor terms. Here, another ten percent might 
be obtained from more extended cluster contributions. 
Third nearest neighbors when added are not converged 
at all. Here a renomalization of at least a factor of two 
is expected from larger clusters. 

More extended clusters and longer range spin correla- 
tions were not covered since alredy the results obtained 
so far led to spin correlation corrections that turn too 
large. This indicates the proximity to the antiferromag- 
netic instability. In some individual cluster calculations, 
already next nearest neighbor contributions turned out 
to be almost as large as nearest neighbor terms. Nev- 
ertheless, for the included clusters the expected compu- 
tation breakdown was not yet seen. Such a breakdown 
must arise as soon as a cluster size is reached where in 
HF-approximation a broken symmetry ground state is 
preferred. Consequently it is concluded that even for the 
half-filled case antiferromagnetic order parameter fluctu- 
ations need to extend over domains that are larger than 
the cluster explicitly included in the computations. 

The longer range correlations contribute differently 
from the nearest neighbor correlations. The direct cou- 
pling matrix elements (O u H) are negligible. However, 
the longer range operators are more sensitive to electrons 
very close to the Fermi surface. Consequently, spin cor- 
relations may form that are dominated by the band elec- 
trons very close to the Fermi energy. 

It is of interest to compare these correlation functions 
with the ones of a two-dimensional Heisenberg model 
p9| . For the latter, the on-site and neighbor correlation 
functions are also given in table IV together with the long 
range limit. It can be seen that the short range corre- 
lation functions of the Heisenberg model are larger than 
the ones of the real system. This is because in the real 
system the Cu3(i :I .2_j / 2-orbitals are more than half-filled, 
and because they are not perfectly correlated. On the 
other hand, the short range correlations of the real sys- 
tem extend already beyound the long range correlation 
pattern of the Heisenberg model. The limiting correla- 
tion function in the magnetically ordered high-T c mate- 
rials is typically lim,,-^ \S(i, i + v)\ = 0.06 |5Q] . This is 
well below the short range correlation functions obtained 
from the present calculation. 







IV. DETERMINATION OF A MODEL 
HAMILTONIAN 

The ab-initio results of the LA provide sufficient infor- 
mation to unequivocally determine a model Hamiltonian. 
Such a condensation of ab-initio results to a model serves 
multiple purposes. One of them is to forward ab-initio 
informations to computations that can no more be per- 
formed on an ab-initio level but only for a model. In the 
following, this applies to the computation of the doping 
dependency of the properties discussed above that were 
calculated for half-filling. Due to the restrictions of the 
Crystal program, such ab-initio calculations would be- 
come costly. Another purpose is that information about 
particular correlations can be represented in form of ef- 
fective interactions. Usually, experiments are fitted by 
models which are represented by the adapted interac- 
tions, but no correlation function for the model is com- 
puted. Such a connection will facilitate comparisons, also 
for differing systems. For the case of the high-T c com- 
pounds, finally, the explicit treatment of correlations was 
so far restricted to models. Therefore, it is of interest to 
see, how well such models match the ab-initio findings. 

The determination of a model from the ab-initio data 
separates into two steps. The first is the choice of the 
model space, i.e. which orbitals to include, and the com- 
putation of the relevant single-particle Hamiltonian. In 
the second step, the effective interactions are computed. 

A. Relevant single-particle space and single- particle 
Hamiltonian 

The information provided by the ab-initio calculation 
concerns atomic orbital degrees of freedom, but so far, 
no more delocalized degrees of freedom like Wannier or- 
bitals extending over sets of atoms. Models built from 
delocalized orbitals can not be directly compared to the 
ab-initio data but would need to be derived in a second 
step from models based on atomic orbitals. 

The smallest model used for the high-T c compounds 
based on atomic degrees of freedom is a 3-band model, 
containing as atomic orbitals the Cu3d x 2_ y 2 and the 
02pb orbitals. However, from the ab-initio charge anal- 
ysis (see table |h|), one notes non negligible fillings of 
Cu4s,p and even deviation from complete filling for the 
02s-orbitals. It is also known that the CuAs orbitals 
contribute actively to the band structure of the most rel- 
evant half-filled band. This was worked out before from 
LDA calculations ]5l| . As a compromise, the selected 
model space for the present application is chosen to con- 
sist of Cu3d x 2_ y 2 and 4s and the 02pi, orbitals. This 
model is also selected because an LDA equivalent exists 

Note that these orbitals can not be seen as a perfect 
representation of the corresponding ab-initio orbitals. 
When taking the charge distribution for the SCF-ground 



state from table |l], then the orbitals included in the 
model represent 4.84 electron per unit cell for the SCF- 
ground state and 4.89 for the correlated ground state in- 
stead of 5.0 as they do for the half- filled band case in the 
model. Consequently, there can be no perfect agreement 
between such a 4-band model and the real system. 

Instead, the most relevant properties are to be 
matched. Here the following properties are selected. The 
first is the half-filling of the uppermost band. This fixes 
the model charge at 5 electrons per unit cell. The sec- 
ond is the exact occupation of the Cu3d x 2_ y 2 orbital ob- 
tained from the respective ab-initio calculation. Since 
the influence of particular correlations will be investi- 
gated, also ab-initio calculations with only partial inclu- 
sion of correlations and varying Cu3c? a; 2_ l) 2-occupation 
will be fitted. The third is the form of the Fermi sur- 
face. The model Fermi surface shall match the ab-initio 
Fermi surface. This is important when longer range cor- 
relations are concerned, and puts restrictive bounds on 
the Cu4s-occupation. A model Cu4s-occupation taken 
from the ab-initio result would lead to a Fermi surface 
that deviates too strongly from nesting. Consequently, 
the Cu4s model occupation is set to 0.25. The omit- 
ted Cw4s-charge in the ab-initio calculation apparently 
stems from bands omitted in the model. Fixing the two 
other occupations freezes the 02p-occpuation. It turns 
out that the deviation of the latter from complete filling 
is only half as large for the model as for the true ground 
state, indicating the bias and the limits of the 4-band 
model. 

Having determined the model occupations for a partic- 
ular ab-initio fit implicitly defines the diagonal or crystal 
field terms e.; of the model. These contain exchange 
and correlation contributions of the omitted degrees of 
freedom as well as exchange contributions due to the 
added on site interactions of the model. Consequently, 
they differ for every fit. Each time, they are determined 
selfconsistently for the model calculation so that the in- 
tended charge distribution is obtained. 

The second set of parameters describes the derealiza- 
tion of the electrons. It consists of the hopping terms. 
Here, it is assumed that the omitted external degrees of 
freedom have no influence on these terms. Since only on 
site interactions will be included, no nonlocal exchange 
contributions arise in the model. The hopping terms 
therefore represent the non diagonal Hartree or alterna- 
tively LDA matrix elements. There are only two relevant 
hopping parameters for this 4-band model. From the 
LDA-fit, it holds for the 3d,2p hopping, tdv — — 1.6ev 
and the 4s, 2p-hopping, t sp — — 2.3eV |plj . At present, 
these hopping elements can not be directly computed by 
the LA since Crystal92 [|J does not separate kinetic en- 
ergy (plus Hartree) terms from exchange terms. How- 
ever, an estimate can be made. The Fock matrix ele- 
ments, fij — tij + Vij(exch), can be computed. For them, 
it is found that fd P — — 2.8eV and f sp = — 4.2eV. When 
approximately correcting for the exchange using the re- 
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lation Vij (exch) ~ — ^ 1 g^ Pij, values of td p = — 1.7eV 

and t sp — — 2.9eV result. The t p d term equals the LDA 
value. This indicates that in the future hopping terms 
between orthogonalized atomic orbitals for models can 
be directly computed by the LA. On the other hand, the 
t sp from the-LDA fit is 20 percent smaller. This differ- 
ence probably arises because the LDA-fit accounts for the 
omission of the 4p-orbitals which otherwise would change 
the higher energy band, while the estimates from the LA 
represent the original bare hopping. In the following, the 
LDA values for the hopping are taken. 

The charge distribution of the original single-particle 
LDA model Hamiltonian (with the LDA values for the e, 
IpH ) amounts to a CM4s-occupation of 0.3, also strongly 
reduced due to the need for an adequate Fermi surface. 
The Cu3c? a; 2_j / 2-occupation for this model is 1.4, and ap- 
parently represents the Cu3d K 2_ y 2-contributions of the 
four included bands. 



B. Interaction terms 

For the determination of interaction terms in the 
model, correlation functions are available. In the model, 
only atomic or on-side interactions are included. These 
are the diagonal interactions for electrons in the same 
orbitals, U3d,U4s,U2 P and the interaction between the 
Cm4s and Cw3(i :I .2_ y 2-orbitals on the same atom, U± S £d- 
These interactions are fitted on a one to one basis with 
the help of the corresponding on site correlation func- 
tions. 

The effective local model interactions Ui are indirectly 
generated from the long range Coulomb interaction that 
prevails in the ab-initio calculation. In this process, dif- 
ferent kinds of rescaling occur. One rescaling process is 
called folding. It is a reduction that is not connected 
to screening. When reducing a single atomic fluctuation 
then not the original atomic interaction is measured but 
the difference between this interaction and the residual 
interaction of the electrons shifted in the process. For an 
almost empty band, the residual interaction is zero while 
in zeroth order, for a half-filled band, it is the neigh- 
bor interaction. A more detailled discussion was made in 
Ref. Jlfj]. Another rescaling is due to screening effects, 
and here two sources exist. The one are the degrees of 
freedom not included in the model, and the other are 
correlations also present in the model but not activated 
for the lack of a longer range model interaction. 

In the following we will add correlations stepwise to the 
ab-initio calculation, and will stepwise interpret the ab- 
initio result in terms of modified model interactions. This 
way, the derivation of a local U from the bare Coulomb 
interaction can be quantitatively understood. The end 
product are local model interactions that represent the 
true ground state. The actual computation procedure for 
the model is as follows. The ground state for the single 
particle part of the model Hamiltonian is easily obtained. 



It is used as an input into the LA program like in the 
ab-initio case the Crystal92 results. The correlation cal- 
culation is then performed by the LA program package, 
but with the interactions reduced to the model interac- 
tions. In principle, for the model, a set of calculations for 
different clusters would be necessary as in the ab-initio 
case before. In reality, the fit was performed by match- 
ing only correlations in one, namely the largest cluster 
to the corresponding ab-initio results. This cluster con- 
sists of five active Cu atoms and four active O atoms. 
Every model single-particle calculation and subsequent 
correlation calculation are embedded into a self consis- 
tent cycle in which the respective single-particle energies 
ej of the model are fixed so that the charge distribution 
of the particular correlated state of the model matches 
the ab-initio counterpart. 

When comparing model and ab-initio correlations, an 
additional constraint needs to be taken into account. It is 
that correlation functions can only be directly compared 
when the respective occupations of the ab-initio calcu- 
lation and the model are identical. A similar constraint 
arises from the folding effect on the value of U. This is 
also strongly occupation dependent. While this poses no 
problem for the Cuid x i_ y i orbital, it involves the other 
two whose charges don't match. Therefore, interaction 
terms for these orbitals are determined in an intermedi- 
ate step in which the model charge balance between the 
Cu4s-orbitals and the 02p-orbitals is shifted so that it 
agrees to the ab-initio result for the momentarily treated 
atom. Fortunately, the different interactions don't influ- 
ence each other much, so that no sizable ambiguity arises 
from this procedure. A test can be made by recomputing 
Usd for the different choices. U^d varies by less then 10 
percent. It is largest for the highest 4s occupation be- 
cause then the 4s, 3d screening explicitly handled by the 
model itself is largest. For the other occupations, larger 
fractions of this screening are mapped by a reduced in- 
teraction parameter U^d- 

Sets of parallel calculations were performed in the fol- 
lowing steps. First, correlations were introduced on a 
single atom only, once for Cu and once for O. In both 
cases atomic charge transfers were alternatively allowed 
or blocked to investigate their effect on U. The corre- 
sponding ab-initio calculation was also restricted to cor- 
relations on single atoms only. This leads to unscreened 
effective interactions U when only single correlation op- 
erators each are included. These interactions U although 
unscreened are folded from the Coulomb interaction and 
do no more represent the original atomic interaction ma- 
trix elements. With charge transfer excluded, the results 
are given in the first line of table |v|. Although folded, the 
U^d is not very far from the atomic interaction matrix el- 
ement which is expected to be U a tom — 25eV. The folded 
interaction U2 P on O on the other hand is quite small. It 
compares to similar interactions obtained before for the 
atomic orbitals on C compounds JTt| . 

The values for U^ s and U^ s ^d are very small. The 
folding effect resulting from neighbor interaction contri- 
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tuitions is expected to be relatively largest for the latter 
terms because the 4s-orbital is most extended. Also it 
should be remembered that the model contains effective 
3d — 4s hopping terms that are 20 percent smaller than 
the original values. For weak correlations it holds that 
correlations scale ~ ~. Correcting for this renomaliza- 
tion would lead to an 25 percent enhancement of the U^g. 
It shall be finally reminded that part of the ab-initio 4s 
occupation apparently stems from other bands and might 
be involved in hopping processes in the ab-initio calcula- 
tion that are not represented in the model. Taking this 
into account should lead to an additional enhancement of 
U^g. For comparison, the value for the effective local in- 
teraction in another system with 4s, 4p electrons, namely 
Ge, is 3.1eV fill , which is not too different from such a 
rescaled value. 

The direct correlation between the 4s and the 3d or- 
bital on the Cu site when added has only a very small 
influence on the value of U%d, and the value of the U^ s 
andL^s^d is not very relevant for U^d- This is unex- 
pected. Similar interactions were important in earlier 
applications when the screening between n and a elec- 
trons in C isomorphs or organic compounds was con- 
cerned [jl(|. In these cases electrons in half-filled bands 
were screened by electrons in wider bands that were also 
half-filled. 

Next, local charge transfers are allowed that arise due 
to correlations. For the single Cu-site, this is a charge 
transfer from the 3d into the 4s orbitals. It leads to 
a strong reduction for U^d (see the second line in table 
|y|). Partly, this reduction originates from the change 
in occupation itself because the folding reduction of the 
original atomic interaction due to longer range Coulomb 
terms is largest for half-filled atomic orbitals and 
the change in Cu3d occupation is toward half-filling. 
The dominant contribution is from the 4s, 3d-screening. 
Such a screening was proposed by Herring long time ago 
f32| |. Astonishingly, it comes into play with the help of 
a charge transfer when starting from the uncorrelated 
ground state, and the major role of the 4s orbitals is 
to represent a reservoir of states. This is interpreted 
in the following way. In HF-approximation, the fluc- 
tuations are somewhat reduced by a too large occupa- 
tion as compared to the one of the true ground state 
(see the discussion of the magneto volume effect in sec- 
tion III C| ). The correct occupation is then obtained by 
the added single-particle correlation operators. This way 
less fluctuations are induced than would arise within a 
uncorrelated ground state that has the correct occupa- 
tion. There are residual fluctuations, part of them from 
the original uncorrelated ground state, and part from the 
charge transfer due to single-particle correlation opera- 
tors. These fluctuations are reduced with the help of 
the two-particle operators. In the model, a sizable part 
of the screening electrons is no more included. When 
fixing the 3d-occupation at the correct value, then very 
much larger fluctuations arise for the model single parti- 
cle ground state than for the original ab-initio uncorre- 



lated ground state. With the original U3d, a suppression 
of fluctuations for the model would result that is very 
much larger than the suppresion of charge fluctuations 
due to two-particle correlation operators in the ab-initio 
case. Consequently, a false description of the 3d correla- 
tion corrections would be made. When the reduction of 
fluctuations is adjusted to the known correction a very 
much reduced effective interaction U^d is found. 

This scenario for the 4s, 3(i-screening is very different 
from the previously treated cases. As just mentioned, 
the screening of the electrons in the half-filled 7r-bands 
in organic systems or C-isomorphs due to the electrons 
in the half-filled a bands was not at all connected with 
a charge transfer but originated solely from two-particle 
correlation operators, i.e. can be seen as a kind of clas- 
sical screening. 

For O, no on-site charge transfer effect did occur. It 
is plausible that the only possible charge transfer which 
is with the not completely filled 2s-orbitals is marginal. 
In a subsequent calculation, all on-site correlations in the 
cluster were treated at once, and also charge transfer be- 
tween the atoms was allowed. The result is displayed in 
the third line of Table [v| There is an additional small 
charge transfer out of the Cu3d-orbitals, leading to a 
small further reduction of U^d- From this term on, the 
U^d obtained for the minimal 4s occupation are given. 
Changing to this reduced value of the 4s occupation con- 
tributes a reduction of U^d of 0.8 eV. This reduction 
arises because the residual 4s occupation is connected 
with a smaller screening. There is also a large charge 
transfer into the 02p orbitals. This causes a reduction 
of the folding effect and a sizable enhancment of U2 P - 
Remember that on this level no other but on-site corre- 
lations were included in the ab-initio treatment. In par- 
ticular the long range part of the Coulomb interaction 
was not screened at all. When nn Cu — O correlations 
are also included (line 4 of Table (vj), both model inter- 
actions are further reduced. The effect is larger for U^d, 
where a charge transfer and screening come together, and 
smaller for U2 P , where a screening gain is reduced by an- 
other enhancment due to the inverse charge transfer. A 
sizable further reduction of U^d occurs when all corrrela- 
tions in the cluster are included. This is mostly due to the 
Cu — Cu spin correlations and the 3d, 2p charge transfer 
caused by them. For U2 P , no value was computed for this 
case. A computation of U2 P requires a strong charging of 
the 4s orbital. In this particular case, it is so large that 
the 3d orbitals are overscreened. This means that even 
for sizably enhanced U^d, too small correlation correc- 
tions of the 3d charge fluctuations were obtained. When 
keeping the original U^d obtained for the other distribu- 
tions, a value of U2 P — 9eV was obtained. This represents 
a further enhancment due to the charge transfer into the 
02p orbitals. 

The values for the model interaction parameters ob- 
tained this way can be seen as upper limits for the true 
model parameters. This is because only specific corre- 
lation corrections were so far included. Part of the ex- 
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eluded correlation corrections can be estimated. This is 
first the effect of the very short range intra-atomic cor- 
relations. Their inclusion led to a reduction of the U by 
10 percent or 0.5 eV for the C-isomorphs |l6j. A sim- 
ilar reduction is expected for O. For Cu a somewhat 
larger reduction is expected due to the higher density 
and due to the screening influence of the filled 3s and 
3p shells. This is next the omitted long range and po- 
larization corrections. Such corrections are irrelevant for 
half-filling. Consequently, they should have no large ef- 
fect on the Usd- However, they should at least reduce 
a large fraction of the enhancement of L^p with charg- 
ing. The resulting estimate of the model parameters is 
given in the fifth line of Table [v|, togehter with a rough 
error estimate. Note that this is the first calculation of 
its kind for an ionic compound as well as for a metal- 
lic compound containing transition-metal atoms. Future 
applications will certainly reduce the uncertainties of the 
error estimates. 

Values for the effective atomic interactions were so 
far computed by LDA frozen charge calculations. From 
these, an effective local interaction U is obtained that 
does not distinguish between different d-orbitals and the 
total angular momentum of the atomic charge. In a sec- 
ond step, higher order Slater parameters were added that 
are taken from experiments on atoms. An introduction 



is given in Ref. |5J 



Table |v| contains an average over 
LDA results [[2l|-|23|, 53) for the resulting diagonal interac- 
tion of the 3d x 2_ y 2 orbitals, U^d, which here is identical 
to the diagonal interactions of the other 3c?-orbitals. Par- 
tially, those calculations also contained results for neigh- 
bor Coulomb interactions V (typically 1 eV or smaller). 
In such cases, the values presented in table [v| are the 
differences U^d — V . The LDA interaction is consider- 
ably larger than the one found from the LA calculation, 
and would result in too large correlations if used for the 
model. 

Note that the definition of U is totally different in the 
two approaches. The one (LDA) freezes charges and does 
not care for their dynamics, i.e. whether they are essen- 
tially localized or whether they are delocalized. Also, 
only the nearest neighbor environment matters. So for 
CuO, almost the same interaction is obtained as for the 
high-T c materials |53]| . The other method (LA) maps all 
particular correlation effects even due to longer range in- 
teractions of the delocalized electrons into an effective 
folded local interaction U. When looking for the deriva- 
tion of the LA value, then it is seen that for the consid- 
ered system, a very peculiar Cu neighbor interaction (or 
spin correlation) leads to a reduction of 2.5eV below the 
LDA- value. Such an effect is 3d x 2_ y 2 specific and would 
not be expected to play a role for the other 3d-orbitals. It 
would also not be expected for CuO for a lack of neighbor 
Cu coupling. For the latter compound, the U^d of the LA 
are expected to be in the range of the LDA values. 

The LDA results were used to explain photoemission 
experiments for different transition- metal compounds. 



Photoemission spectra calculated for an Anderson im- 
purity model with the computed LDA values of [73(2=9.5 
eV pa ] led to very good agreement with experiment for 
CuO and (with a particular exception) for Nd,2CuO^. 
A model when mapped to the same experiment led to a 
value of C/ 3d =8.4 eV for CuO @- This demonstrates 
that the LDA values are the correct values for the inter- 
actions among the completely localized 3d-orbitals. How- 
ever, this can not be taken as evidence for the correctness 
of the value for the 3(4,2 _ y 2 orbitals. A modification of 
the latter interaction towards the LA result would prob- 
ably not change the computed spectra very much. There 
is a small deficiency for Nd2CuO±, though, when fitted 
to the LDA values. This is the existence of a local sin- 
glet peak at the upper band edge in the calculation, which 
also shows up in calculation and experiment for CuO, but 
not in the experiment for Nd2CuO^. It is a valid specula- 
tion whether already the reduced 3d x 2_ y 2 U^d originating 
from the effective neighbor Cu-correlations in the plane 
would remove this deficiency. 

The U2 P interaction of the LA is in qualitative agree- 
ment with earlier values found for the 2s, 2p-interactions 
in diamond (7.2 eV [^l)) but larger than the LDA esti- 
mate. There exists a spectroscopical fit for L2 P =5.5eV 
@. 

The presented difference between the LA values and 
the LDA values of U^d matches the difference between 
LDA and experiment found earlier for the transition met- 
als. For the transition metals, the U^d of the LDA are 
apparently independent of band filling 56 57]] while the 
experimentally needed quantities are strongly filling de- 
pendent and considerably smaller - except for the com- 
pletely filled 3d-band limit |IJ,|l^,|l^,|l9| . It was proposed 
to resolve this discrepancy and the filling dependence of 
U^d by a not fully screened neighbor interaction V. This 
would explain the LDA deviations and the filling depen- 
dency of U3d as a folding effect fl6|] . 

Calculations with the LA are now feasible for transition 
metals. The results obtained here give hope that from 
such calulcations, apropriate values for the U^d of the 
transition metals can be obtained. For comparison, the 
value of U3d needed for Ni is 4.7eV, and not very much 
lower than the final LA estimate for Cu in SrCu02- 



V. MODEL INTERACTION AND SPIN 
CORRELATIONS 

The model interactions derived in the last section were 
optimized with respect to charge distribution and on-site 
correlations. Next, we will control how well this model 
is also able to reproduce the most interesting longer 
range correlation features, namely the spin correlations 
between different Cu3c?-orbitals. 

For the 5 Cu-cluster, the results of the model calcula- 
tion with U^d — 6.3eV are compared to the corresponding 



13 



ab-initio values in table VI. As can be seen, the model 



neighbor spin correlations are only half of the ab-initio 
values, and the 2nd neighbor terms are reduced to one 
sixth. Consequently, the 4-band model with only on-site 
interactions included can not consistently reproduce on- 
site correlations and neighbor spin correlations. Ab-initio 
neighbor spin correlation functions are only matched by 
a U^d that is 50 percent to large, and even then, the defi- 
ciency of the second neighbor correlation functions is not 
completely removed (see table VI). 

Part of these discrepancies can be understood by con- 
sidering the U dependence of the individual correlation 
functions calculated by means of the weak-correlation ex- 
pansion. Here, the proximity to the magnetic instability 
plays a particular role. For the model (and for the ab- 
initio calculation), the interaction enters in two quanti- 
ties, namely {OH) and (OHO) c . If the (OHO) c were not 
depending on U , then the variational parameters (and 
the correlation functions up to saturation) would rise lin- 
early with interaction strength. The model results are 
different and indicate that the interaction dependency 
of the second terms must come into play. This holds 
particularly true for the longer range correlations who 
rise very much more than linearly with interaction. This 
anomalous U dependency can only be understood by the 
proximity of the magnetic phase. As discussed before, 
the terms (OHO) c represent the two particle excitation 
energies. Close to a magnetic instability, these might be- 
come very much smaller and tend to zero, leading to an 
anomalous U dependence of the correlation parameters. 
Apparently, this applies to the model. 

This interpretation also explains why the longer range 
spin correlations in the model are relatively weaker than 
for the ab-initio case. The diagonal terms of (OHO) c 
represent energy differences of bare excitations out of 
the SCF-ground state wave function. This means that in 
the ab-initio calculation, for these matrix elements the 
uncorrelated or HF susceptibility enters. In the model 
however, the calculations were not performed with bare 
but with screened interaction parameters. This means 
that for the longer range spin correlations the energy 
difference of bare excitations is computed with screened 
interactions and therefore contains correlations to some 
extend, in contrast to the ab-initio calculation. Conse- 
quently, the model result for the longer range spin cor- 
relations might be more adequate than the ab-initio re- 
sult and might even indicate by which amount the ab- 
initio results need to be corrected. With only second or 
third neighbor correlations included, the model is still 
far from instability, in contrast to the conclusion derived 
from the ab-initio calculation. A magnetic instability 
might only occur when considerably longer-range mag- 
netic correlations are added. It can even not be excluded 
that the charge transfer connected with magnetic corre- 
lations causes a first order phase transition, and that no 
divergency of the long-range correlations can be detected 
in the metastable state without broken symmetry. 

Nearest-neighbor spin correlations don't display such a 



strong {/-dependence. While it can not be excluded that 
the difference between ab-initio and model results might 
also originate from the overestimated Stoner enhance- 
ment, there is another deficiency of the model that points 
to a different source. The neighbor spin correlations are 
connected to an explicit Cu3d-02p charge transfer that 
is by almost one order of magnitude smaller in the model 
than in the ab-initio calculation. A reason for this may 
be that the CuAs,p degrees of freedom are mostly re- 
moved from the model and only indirectly included in 
the form of reduced on-site interactions. It might well be 
that these omitted degrees of freedom contribute more 
actively to the neighbor-Cu-spin correlations with the 
help of an induced magnetic exchange interaction be- 
tween neighbor Cu sites. Another reason for the discrep- 
ancy between ab-initio and model results might be that 
the O occupation in the model is considerably larger than 
in the ab-initio case, and might also reduce the charge 
transfer. 

Apparently, the 4-band model with on-site interactions 
only is not quite adequate to deal with the most inter- 
esting outcome of the ab-initio calculations, namely the 
anormalous neighbor spin correlations. 

The Stoner enhancement in the longer-range spin cor- 
relation functions of the model calculation of the 4 (or 
3) band model and also apparently in the ab-initio cal- 
culation is very different from results expected for a 
single band Hubbard model. When using a weak cor- 
relation expansion, then it is well known that for the 
one-dimensional case all interaction contributions in the 
terms (OHO) c drop out. This means a linear rise of 
all correlations with U but no Stoner enhancement. For 
the two-dimensional model with perfect nesting, simi- 
lar results are expected. It can not be ruled out that 
closed loop terms lead to interaction contributions in the 
(OHO) c for the non nesting case but such terms are not 
yet present in the considered 5 atom cluster. This in- 
dicates that the magnetism in the real system is essen- 
tially of itinerant or spin-density- wave nature ( although 
strongly enhanced by the almost perfect nesting), and 
that a simple single-band Hubbard model might not to 
be the correct approximate description. 



VI. THE MODEL AWAY FROM HALF-FILLING 

The computations performed so far were restricted to 
the so called half-filled band case. As mentioned before, 
the program Crystal92 can only be used for an integer 
number of electrons per unit cell. There is no such re- 
striction for the LA program package. Consequently, the 
model calculations can easily be extended to partial fill- 
ings. For simplicity, the model SCF calculations were 
not repeated for differing fillings but the single-particle 
Hamiltonian at half-filling was frozen in, and only the 
Fermi energy was shifted. This approximate treatment 
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seems justified because contributions relevant for charge 
redistribution like the long-range Madelung terms are not 
included in the simple on-site interaction model. 

Of interest is the dependency of the neighbor and of 
the longer range spin correlations on band filling. Fig. |^ 
displays the nearest-neighbor [y = 1) and next-nearest- 
neighbor (y — 2) Cu-spin correlations as a function of 
the occupation of the uppermost band {jib)- These cor- 
relation functions were taken from calculations for a five 
Cu atom cluster again. Corections towards the full re- 
sults are typically 20 percent for the nearest-neighbor 
terms and more than 100 percent for the second-nearest- 
neighbor contributions. Two values for the interaction 
parameter U^d were taken, namely the value deduced 
from the ab initio fit (6.3eV, continuous lines), and a 
value enhanced by 20 percent (7.8eV, dotted line). The 
second computation with an enlarged U%d was made to 
obtain an estimate for the 4-band-model shortcomings in 
comparison to the fictituous ab-initio result. 

Both correlation functions reduce in strength when 
electrons are removed. The longer range function does 
so somewhat stronger. However for the range of interest, 
i.e. around optimal doping (rid — 0.8), both functions 
are still sizable and not very much smaller than for the 
metastable non magnetic half-filled case. This a posteri- 
ori justifies the choice of such a metastable state in the 
ab-initio correlation calculations. It also demonstrates 
that for all fillings of interest very sizable nearest neigh- 
bor short range antiferromagnetic correlations exist to- 
gether with longer range itinerant antiferromagnetic fluc- 
tuations. The {/-dependence of the longer-range correla- 
tion function is more pronounced and is strongest close 
to half-filling, indicating again the underlying Stoner en- 
hancement. 

Also of interest is the change of the charge distribu- 
tion with refilling. In the single-particle approximation, 
the electrons close to the Fermi surface are mostly d-likc. 
Fig. H displays the non d-fraction of the density of states 
as a function of doping. It can be seen that close to the 
Fermi energy it amounts to 10 percent. As mentioned be- 
fore, the model itself is unable to account for Madelung 
corrections that would certainly modify such an extreme 
density distribution of the removed charge. Also a possi- 
ble redistribution that would come in with the self con- 
sistent computation is not included. However, we will 
discuss investigated to which extent correlations lead to 
a redistribution of the removed charge. As can be seen 
from Fig. [| there is indeed a sizable change of charcter 
of the removed density. Most of this is change into 02p- 
character. Again, the computation is performed for two 
values of U-^d- The smaller value (6.3eV) leads to a 20 
percent charge redistribution while the larger value leads 
up to 50 percent corrections close to half-filling. Photoe- 
mission experiments not too far from the magnetic state 
found indeed that the electrons removed from the system 
were largely of 02p-character |5£| . 



VII. SPIN CORRELATIONS AND NEUTRON 
SCATTERING RESULTS 

In the ab-initio calculation for the half-filled case, very 
strong antiferromagnetic neighbor Cu spin correlations 
in connection with a Cu — O charge transfer and with 
longe range antiferromagnetic polarizations were found. 
The subsequent model calculations have shown that the 
neighbor spin correlations are not restricted to the im- 
mediate vicinity of half-filling but exist for every filling. 

This prediction can be tested by comparing the cal- 
culated results with quantitative magnetic neutron mea- 
surements. From these experiments, a quasi-equal time 
spin correlation function S(Q) was obtained for the 
metallic compound LaQ.^SrQ ^Cu20^ by extending the 
energy integration up to 0.45eV. The data show a strong 
longer range structure that is expected to exist indepen- 
dently of an also found small quasi-elastic scattering aris- 
ing from an incommensurate spin density wave present in 
this particular compound. In the following, a comparison 
is made between the theoretical equal time correlation 
function and the measured quantity. Both are not iden- 
tical. The theoretical quantity is obtained for the infinite 
layer case where for 0.15 holes per unit cell apparently no 
incommensurate spin wave exists, and consequently less 
magnetic scattering is expected than for the measured 
compound. On the other hand, the theoretical quantity 
represents the true equal time correlation function and 
contains contributions that are not in the range of the 
measurement. 

The limiting equal time case of the measured correla- 
tion function is defined as 

S(Q) = -L J d 3 r J dVS(r, ^e^""' . (21) 

The theoretical spin correlation function is derived from 
the model calculation at apropriate doping. It is repre- 
sented by spin correlations between different orthogonal 
orbitals, 

S(i,j,G) = (V corr \s i (0)s j (G)\* corr ) . (22) 

Here i denotes the i-th orbital in the unit cell with atom 
position fi, G represents the lattice vectors, and sl(G) 
represents the spin operator for orbital i in the G unit 
cell. When assuming that the spatial moment distribu- 
tion is shrunk to the nuclear positions, 

S(r,r*)~ ]T d(r-n + G)5(f -? j + & + &)S(i,j,&) 

(23) 

one obtains 

S(Q)=J2S(i,j,d)e i ^- r ^ . (24) 
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This function is very easy to compute. Fig. Q contains 
the results for a particular Q direction, namely the di- 
agonal (1,1) axis, obtained in different approximations. 
The zone boundary is at h=l, the intensity is given per 
formular unit which here is equivalent to a unit cell or to 
a single Cu atom. 

The lowest curve represents the result for the single- 
particle ground state. It represents the exchange holes. 
As the finite value at h = indicates, the summation 
in eq. [24| was not brought to convergency. Instead, the 
Cu — Cu density matrix elements were only included up 
to the 4th neighbor, and no density matrix elements with 
4s or 2p orbitals extending beyound the nearest neighbor 
Cu — O terms were added. The maximal deviation occurs 
for h — were the contributions from all missing terms 
add up. Due to dephasing, the correction is very much 
smaller for finite h. Due to the fine structure in the 
unit cell, this function is finite at the first lattice vector 
(h = 2). This represents the Cu — O correlation function. 

Next, short range correlations as they are deduced 
from a single coherent 5 Cu cluster calculation are 
included (dotted curve). Here, the nearest neighbor 
Cu — Cu correlations come into play and cause a peak 
at the zone boundary (h — 1). When extending the cor- 
relation treatment to a 9 Cu cluster, second and third 
neighbor correlations are more correctly treated. They 
lead to a narrowing of the peak and to a small enhance- 
ment (continuous curve). Finally, also the corresponding 
values with enlarged U (7.8eV instead of 6.3eV) are given 
(broken curve). Increasing U leads to a strong enhance- 
ment of the maximum, indicating again the proximity to 
a magnetic phase transition. These results are compared 
to experiment (dots) in Fig. |]. As expected, the the- 
oretical equal time correlation function is always larger 
than the experimental correlation function whose energy 
integration extends only to 0.45 eV. Beyound h=1.5, the 
experimental results are influenced by the next Bragg 
peak, and no more meaningful. 

There are specific contributions to the theoretical cor- 
relation function that are not expected to be seen by 
experiment. These are the short range contributions con- 
nected with the Cu — O hopping, arising already with- 
out correlations. The hopping energy connected with 
this part of the correlation function is t = 1.6eV, and 
very much larger than the energy cut-off. Consequently, 
only a marginal part of these contributions is expected 
to show up in experiment. A considerably larger frac- 
tion of the uncorrelated longer range Cu — Cu contribu- 
tions is expected to show up since these mostly arise from 
the uppermost band. Also the relevant correlation con- 
tributions are expected to be measured by experiment. 
While the on-site correlation functions might not fully 
show up, the effect of the neighbor Cu — Cu spin cor- 
relations is expected to arise mostly from the electrons 
in the uppermost band, and the longer range enhanced 
spin correlations are certainly connected with electrons 
close to the Fermi surface, as is indicated by their strong 



resonance dependence on U. Fig. |^ displays the cor- 
relation functions of Fig. [|, but with all Cu — O con- 
tributions of the single-particle approximation removed, 
and with the residual function shifted so that 5(0) = 
holds. The residual Cu — Cu single-particle contributions 
are small and essentially bell shaped (lowest curve). The 
correlation contributions lead to a pronounced maximum 
around the zone boundary. When correlations result- 
ing from the 9 Cu-cluster calculation are included, then 
the half width of the correlation peak corresponds well 
to the half width of the experimental peak. However, 
for the value of U taken from the fit to on-site corre- 
lations, the integrated scattering intensity is not larger 
than the experimental counterpart. The result for a 20 
percent enhanced U finally leads to a correlation func- 
tion that is systematically larger than the experimental 
curve. The ab-initio calculation if performed for the rel- 
evant doping would certainly give a correlation function 
as large or even somewhat larger than the model result 
for the enhanced interaction. A future comparison with 
experimental results for a metallic compound without a 
spin density wave will allow to decide whether the model 
results or the ab-initio results are more reliable. The 
ab-initio calculations might overestimate the Stoner en- 
hancement, while the model might well leave out relevant 
degrees of freedom, and might consequently need to be 
extended. 

The theoretical results represent not only the partic- 
ular doping of 0.15 holes but should be representative 
for a wider range of doping even farther away from the 
magnetic case. As Fig. || demonstrates, the neighbor 
correlations that represent the weight of the peak around 
the zone boundary reduce only slowly with further dop- 
ing. The longer range correlations are expected to re- 
duce faster, so that a continuous widening of the peak 
with further doping is expected. These correlation fea- 
tures appear over a rather wide range of doping and have 
consequently, no direct connection with any kind of Mott 
Hubbard transition. 

The equal-time spin correlation function was computed 
earlier for a one-band Hubbard model 59 6(]] or for a t- J 
model |6^]. In these computations, due to the particular 
Fermi surface, at 0.15 doping a spin-density wave shows 
up. The resulting equal time correlation function is dif- 
ferent from the one given in Fig. ||. It is close to the 
bell shaped curve of the uncorrelated electrons in Fig. |5| 
but enhanced by a factor of 3. In addition, for magneti- 
cally ordered states, there is a very narrow peak just at 
h = 1.0 ( or a set of two peaks close to this point). This 
peak dissapears for the not ordered states, but its width 
is usually not resolved due the finite k-point mesh used 
in these computations. There is no evidence for strong 
shorter- or longer-range spin-correlation features in the 
nonmagnetic metallic state. This indicates that single- 
band models with local interactions don't adequately de- 
scribe the low-energy degrees of freedom of the metallic 
case. 

The extended range of longer-range antifcrromagnctic 
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correlations found in the ab-initio calculation but also 
for the 4-band model, contrasts to single-band model re- 
sults. This is connected to the following difference . On- 
site correlations are strongest for a 3d x 2_ y 2 occupation 
of 1. This occupation occurs at 0.4 to 0.5 doping. An- 
tiferromagnetic order on the other hand is strongest for 
perfect nesting in the half-filled band case. Inbetween 
both points, a region of strong fluctuations is found. For 
a single-band model, these different points are reduced 
to a single point, half-filling. 



VIII. CONCLUSIONS 

The results obtained within the framework of the LA 
can be hierachically classified into the following cate- 
gories. The first concerns the general field of ab-initio 
calculations for the transition metals, the next deals with 
the connection of full Hamiltonian and model Hamiltoni- 
ans, and the last one finally covers the specific electronic 
properties of the high-T c compounds. 

Concerning the general field of ab-initio calculations 
for transition metals, we have presented the first LA com- 
putation which also provides detailed correlation func- 
tions, and which, as mentioned is not connected to a 
homogeneous-electron-gas like approximation. We could 
explicitly investigate those correlation effects that are out 
of the reach of the LDA. An example relevant for the gen- 
eral field is the correlation induced 3c? — 4s charge trans- 
fer. A similar charge transfer is expected from magnetic 
order in the transition metals and was proposed long ago 
by Lang and Ehrenreich as an explanation for the in- 
verse magneto- volume effect in Ni [37]]. This can now be 
quantitatively adressed. 

We have found that the weak-correlation expansion in 
which the LA is computed can be successfully applied 
to systems as strongly correlated as the high-T c mate- 
rials. Astonishingly, problems arose neither in the con- 
text of strong atomic correlations nor due to a possi- 
ble Mott-Hubbard transition on either the atomic or a 
more extended unit-cell scale, but only from the close- 
ness of the ground state to a magnetic phase and from 
the resulting Stoner enhancement. A future extension of 
this weak-correlation expansion, from the linearized to 
the full CCSD equations in the restricted operator space 
should help to improve this specific shortcoming. 

The calculations have also demonstrated how impor- 
tant it is to use an SCF-calculation for the full solid as 
a starting point. In the past, it has been necessary to 
restrict oneself to calculations for small Cu ~ O clus- 
ters when applications based on ab-initio treatments be- 
yound a homogeneous-electron-gas-based approximation 
were made. One of the first such calculations has also 
been done in the framework of the LA ]34| , and a com- 
parison of the results clearly shows how adversely clus- 
ter constraints influence basic results like the electronic 
charge distributions or correlation functions. 



Concerning the transition from the ab-initio calcula- 
tion to a model Hamiltonian, the dominant issue is the 
determination and the analysis of the interaction pa- 
rameters of the effective Hubbard models. Among oth- 
ers, we have provided a detailed derivation of the effec- 
tive CuZd x 2 _ y 2 -interaction parameters, starting from the 
bare Coulomb interaction, and analysed in particular the 
screening effect of the 4s, 4p-electrons that had been pro- 
posed long ago by Herring [|52| . A surprising finding was 
that this screening is not much connected with the resid- 
ual interactions between the 3d and the 4s, 4p electrons, 
but largely mediated by a charge transfer from the 3d 
orbitals into the 4s, 4p orbitals, when starting from the 
SCF-ground state. 

The obtained Cuid x 2_ y 2 -interaction parameter turned 
out to be somewhat smaller than the global Cu3d- 
interaction parameters that were determined by frozen 
charge LDA calculations. The difference apparently 
results from residual interactions between electrons in 
3d x 2 _ y 2 -orbitals on neighbor-Cu sites that are only ac- 
counted for when correlation functions are used as a 
means to determine the interaction parameter. A sim- 
ilar deviation was noted earlier for the case of the 
transition-metal interaction parameters, for which some- 
what smaller values were obtained from fits to experiment 
than from computations by the LDA Jl(| . 

Finally, we will adress the specific properties of the 
metallic CuO compounds. On the single-electron level, 
our results are similar to the LDA results. This con- 
cerns in particular the relevance of the 4s orbitals for 
the dispersion of the half-filled band and for the form of 
the Fermi surface. However, there arc also differences. 
Surprising is the one for the 3d-occupation which comes 
out too large in the LDA. Furthermore, we found siz- 
able correlations on different length scales. While the 
strong atomic correlations were expected for these com- 
pounds, we found in addition a strong magnetic nearest- 
neighbor Cu — Cu correlation that might even lead to a 
neighbor attraction of electrons with different spin. This 
correlation is not due the longer-range magnetic fluctu- 
ations, however it may well be enhanced by it. It is 
also connected with a sizable Cu — O-charge transfer. 
A homogeneous-electron-gas-based method like the LDA 
is neither expected to be able to handle such a corre- 
lation nor the connected charge transfer. This explains 
the just-mentioned difference in the 3d-occupation. In 
addition, a sizable long-range magnetic polarization was 
found that can best be described in terms of a Stoner 
enhancement. All these features turn out to be present 
over a large doping range, and not only very close to 
half-filling. 

The connection between variations in the magnetic cor- 
relations and the charge transfer is expected to result in 
interesting couplings between the magnetic and lattice 
degrees of freedom. In particular, it should not be sur- 
prising if the magneto-volume effect for the case of the 
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half-filled-band systems turns out to be very small or 
even negative, as is the case in Ni. 

The just-mentioned particular neighbor spin correla- 
tions, enhanced by longer range magnetic fluctuations, 
dominate the spin correlation function, and explain the 
features found in the measured spin correlations for 
£ao.85<SVu 15C112O4 Q. Apparently, neither the single- 
band-Hubbard model results nor the t-J-model results 
can explain th is sp in correlation function (see discus- 



sion in section VII). It was not even possible to bring 
the results of the 4-band model with only on-site inter- 
actions to good agreeement with the ab-initio results, 
when on-site and longer-range correlations were jointly 
concerned. It seems that a proper description of the real 
system can only be obtained if in such a 4-band model the 
background-induced magnetic Cu neighbor interactions 
are explicitly taken into account, or if the model is gener- 
alized by an explicit inclusion of the Ap orbitals, may be 
even of the full screening process of the 3d-interactions 
by the 4s, 4p-orbitals. In our future work, we shall inves- 
tigate such extensions. 

To conclude, ab-initio correlation calculations can now 
be performed for the transition metals. With the local 
ansatz, details of the correlation functions as well as a 
good understanding of the relevant short-range correla- 
tion features can be obtained. The first application for 
a metallic high-T c compound shows a fairly good agree- 
ment between the computed and the measured magnetic 
correlation functions. 
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Correlations 






i? corr (a.u./u.c.) 


spin contrib. 




on site 






-.1248 








ou-u nn 






-.Ulo4 




0.20 




Cu-Cu nn 






-.0190 




0.80 




O-O nn 






-.0033 








Cu-0 nnn 






-.0036 








Cu-Cu nnn 






-.0157 




0.95 




Cu-Cu nnnn 






-.0040 




0.95 




TABLE I. Correlation energy contributions in 


atomic units for particular 


succesively added operators. 


Orbital 


HF 


on site corr 


nn corr 


full corr 


LDA 


Cu3d x 2_ y 2 




1.51 


1.33 


1.17 


1.15 




Cu3d z 2 




1.95 


1.94 


1.94 


1.94 


9.30* 


Cu3d X y , 3d xz 




2.00 


2.00 


2.00 


2.00 




CuAs 




0.50 


0.55 


U.O I 


0.58 


0.53 


CuAppi 




0.30 


0.33 


0.34 


0.34 


0.64* 


CuAp± 




0.09 


0.11 


0.11 


0.11 




02s 




1.82 


1.82 


1.81 


1.81 




02p b 




1.42 


1.48 


1.57 


1.58 




02p ort h 




1.97 


1.96 


1.96 


1.96 




02 P± 




1.95 


1.93 


1.92 


1.91 





TABLE II. Charge distributions for the SCF ground state and with correlations added, in comparison to LDA results [|42| 
* The sums over five respectively three partial contributions each are given. 



Orbital 


rii 


A(HF) 


A(corr) 


A ma x(corr) 


Cu3d x 2_ y 2 


1.17 


0.243 


-0.122 


-0.172 


CuAs 


0.57 


0.203 


-0.009 


-0.081 


CuAp p i 


0.34 


0.141 


-0.005 


-0.029 


02s 


1.81 


0.086 


-0.002 


-0.009 


02p b 


1.57 


0.168 


-0.033 


-0.049 



TABLE III. On-site correlations for the different atomic orbitals. The individual terms are defined in the text. 



Included correlations 


v = 


v = 1 


v = 2 


v = 3 


HF ground state 


0.276 


-0.012 


0.001 


0.001 


on site 


0.478 








up to v = 1 


0.530 


-0.140 






up to v = 2 


0.540 


-0.220 


0.170 




up to v = 3 


0.543 


-0.243 


0.183 


0.069 


Heisenberg model 


0.75 


-0.34 


>0.10 


>0.10 



TABLE IV. Spin correlation functions for C , M3d a .2_ !/ 2-orbitals between neighbor sites i,i + v, as functions of the included 
correlation operators, in comparison to the 2d-Heisenberg model. 
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Correlations 


u d 




Us 


u sd 


u p 


on site, single without charge transfer 


20.8 




1.8 


1.3 


D.D 




nn ci + p cincrlp 


12.0 




1.7 


1.1 


6.6 




on site, global 


10.4 




1.8 


1.2 


8.0 




CuO-nn, global 


8.8 




1.8 


1.1 


7.6 




all, global 


6.3 




1.8 


1.1 






estimate 


5.7 


±1.0 


1.8 


1.1 


6.0 


±1.0 


LDA 


9.0 


±1.0 






4.5 


±2.0 



TABLE V. Effective on-site interaction parameters Ui (eV), as obtained in different approximations. The final estimate is 
also given in comparison to typical LDA results fell B3U53] . The individual terms are explained in the text. 



Computation 


A 3d 












v = 1 


v = 2 


v = 3 


ab initio 


-0.122 


-0.140 


0.075 


0.072 


U d = 6.3eV 


-0.122 


-0.071 


0.013 


0.011 


U d = 7.8eV 


-0.145 


-0.097 


0.025 


0.024 


U d = 9.4eV 


-0.160 


-0.126 


0.045 


0.043 



TABLE VI. On site and ^th neighbor Cu — 3d x 2_ y 2 correlation functions for the cluster with 5 active Cw-atoms. Ab initio 
results in comparison to model results with differing U d 



FIG. 1. Schematic representation of the two largest clusters for which interaction matrix elements for the basis orbitals were 
generated. Cu-atoms are denoted by crosses, O-atoms for which correlation operators were included are denoted by filled 
circles, while O atoms who contributed only to the Vijki are denoted by open circles. 

FIG. 2. Nearest neighbor (u — 1) and second nearest neighbor (y = 2) Spin correlation function in dependence of the filling 
n b of the uppermost band, obtained for a model with U3 d =6.3 eV (continuous lines) or 7.8 eV (broken curve), respectively. 

FIG. 3. Relative non 3d-like density of states at the Fermi energy in dependence of the filling rij, of the uppermost band, 
obtained without correlations (dotted-broken curve), for U3 d =6.3 eV (continuous lines) and 7.8 eV (broken curve), respectively. 

FIG. 4. Equal time spin correlation function S(Q) for Q = (h,h,0) in comparison to experiment (empty circles). Given 
are the results of the HF-ground state (broken-dotted curve), the 5 atom cluster result (dotted line) and the 9 atom cluster 
result (continuous line) for [/3d=6.3eV, and the 9 atom cluster result for L?3d=7.8eV. 



FIG. 5. Cu — Cu -dependent part of the equal time spin correlation function S(Q) in comparison to experiment. Definition 
of the curves like in fig. H 
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